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Abstract 

We consider a possible generalization of the random matrix theory, 
which involves the maximization of Tsallis' q-parametrized entropy. We 
discuss the dependence of the spacing distribution on q using a non- 
extensive generalization of Wigner's surmises for ensembles belonging to 
the orthogonal, unitary and symplectic symmetry universal classes. 

PACS numbers: 03.65.-w, 05.45. Mt, 05.30.Ch 

1 Introduction 

Non-extensive statistics is by now recognized as a new paradigm for statis- 
tical mechanical consideration yQ. It revolves around the concept of Tsal- 
lis's information-entropy measure S qi a generalization of the Boltzmann-Gibbs- 
Shannon (BGS) entropy, that depends on a real parameter q. This parameter, 
which is called entropic index, characterizes the degree of extensivity of the sys- 
tem. If the system is composed of two independent systems A and B such that 
the probability p(A + B) = p(A)p(B), the entropy of the total system is given 
by the following pseudo-extensive relation 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B), (1) 

from which the denunciation non-extensive comes. Tsallis 's entropy becomes 
Shannon's measure for the particular value q = 1. BGS entropy is the natural 
entropy for modelling chaotic dynamics that suffer from exponential sensitivity 
to initial condition, having at least one positive Lyapunov exponent. Tsallis 
statistics provide possible way to describe non-integrable systems when the sen- 
sitivity to the initial conditions diverges less than exponentially. This has been 
demonstrated by studying a great variety of systems. Among them is the logistic 
map 01 El EH, which have been used since more than two decades to demonstrate 
various routes to chaos 0. 

Quantum chaotic systems are known to be described in terms of the random- 
matrix theory (RMT). It models a chaotic system by an ensemble of random 



Hamiltonian matrices H that belong to one of the three universal classes, or- 
thogonal, unitary and symplectic. The Gaussian probability density distribution 
of the matrix elements P (H) oc exp [— rfli [H< H jl is obtained by maximizing 
BGS entropy [21 EH ■ A complete discussion of the level correlations even for 
these three canonical ensemble is a difficult task. Most of the interesting re- 
sults are obtained for the limit of N — > oo. Analytical results have long ago 
been obtained for the case of TV = 2 |3| . It yields simple analytical expressions 
for the nearest-neighbor-spacing (NNS) P(s), renormalized to make the mean 
spacing equal one. The spacing distribution for the Gaussian orthogonal en- 
semble (GOE), P(s) — -|se~~ s , is known as Wigner's surmise. It is obtained 
by imposing the requirement of being invariant with respect to the base trans- 
formation 8,9. Analogous expression are obtained for the Gaussian unitary 
and symplectic ensembles (GUE and GSE), respectively [HI EL an d given here 
below. These expressions have been successfully applied to the NNS of various 
chaotic systems. On the other hand, there are elaborate theoretical arguments 
by Berry and Tabor ^U] that NNS of classically integrable systems should have 
a Poissonian statistics. The Poisson distribution exp(— s) of the regular spectra 
has been proved in some cases (see, results by Sinai an d Marklof ^5] > f° r 
instance). Still its mechanism is not completely understood. It has also been 
confirmed by many numerical studies, e.g. |13| . so that the appearance of the 
Poisson distribution is now admitted as a universal phenomenon in generic inte- 
grable quantum systems. Nevertheless, it is well known that not all the regular 
systems have a Poissonian NNS distribution. The two-dimensional harmonic 
oscillator is a classical example ^01- In this case, numerical calculations show 
the NNS distribution vanishes everywhere except in a narrow region centered 
at s w 1. This suggests that the spectrum has nearly a picket-fence shape. 

In this work, we suggest that the non-extensive generalization to RMT pro- 
vides a principled way to accommodate systems with mixed regular-chaotic 
dynamics. We use the maximum entropy principle, with Tsallis' entropy, to 
obtain the distribution functions for random-matrix ensembles belonging to the 
three canonical symmetry universalities. 



2 Nonextensive RMT 

The Tsallis entropy is defined for the continuous variables entering the Hamil- 
tonian matrix H by 

where 

0-1 

dH=n <m£ii[ n dH-w, (3) 

m>n 7—1 7n>n 

while = 1,2 and 4 for the Hamiltonians having orthogonal, unitary or sym- 
plectic symmetries. We shall refer to the corresponding ensembles as the Tsallis 
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orthogonal ensemble (TsOE), the Tsallis Unitary ensemble (TsUE), and the 
Tsallis symplectic ensemble (TsSE). For q — > 1, S q tends to BGS entropy 
Si [Pp{q,H)\ = — J dHPp(q,H)lnPp(q, H), which yields the Gaussian ensem- 
bles QIIHI- There are more than one formulation of non-extensive statistics which 
mainly differ in the definition of the averaging. Some of them are discussed in 
|14| . We apply the most recent formulation |T2]- The probability distribution 
Pfj{q, H) is obtained by maximizing the entropy under two conditions, 

dHP p (q,H) = 1, (4) 
= < T ~ (5) 



fdH[Pp(g,H)] q Tr(HiH) 2 
jdH[P (q,H)]« a ? 



where erg is a constant. The optimization of S q with these constraints yields a 
power-law type for Pp(q,H) 

Pp(q,H) = Z- 1 [1 + (q - l)\ {Tr (H+ff) - 4}]'^ , (6) 

where r\ q > is related to the Lagrange multiplier 77 associated with the con- 
straint in (5) by 

% = V/ J ' dH[Pp{q,H)] q , (7) 

and 

Z q = JdH[l + (q- 1)17, {Tr (rfH) - <J 2 p }] . (8) 

It turns out that the distribution (6) can be written hiding the presence of er^ 
in a more convenient form 



where 



and 



P p (q,H) = Z- 1 [l + (q-l) Vq TT(H^H)y^ , (9) 

"« = /dff[P^(g,ff)]« + (1-^' ^ 

^ = JdH[l + (q- l) Vq Tr (fftff)]-^ _ (n) 

The non-extensive distribution (9) reduce to the statistical weight of the Gaus- 
sian ensemble when q = 1. 

We now calculate the joint probability density for the eigenvalues of the 
Hamiltonian H . With H = U^XU, where U is the global unitary group. For 
this purpose, we introduce the elements of the diagonal matrix of eigenvalues 
X = diag(xi, • • • , xjv) of the eigenvalues and the independent elements of U as 
new variables. Then the volume element (3) has the form 

dH = \A N (X)fdXdn(U), (12) 
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where A a? (X) = Y[ n>m (x n — x m ) is the Vandermonde determinant and dfi(U) 
the invariant Haar measure of the unitary group In terms of the new vari- 
ables, Ti[H'H) = Yli=i x i so ^ na t the right-hand side of Eq. (9) is independent 
of the "angular" variables in U. Integrating (6) over fi(U) yields the joint prob- 
ability density of eigenvalues in the form 



,Xn) 



n ( Xr > 







1-(1-?)t? 9 



N 



(13) 



where Cp is a normalization constant. This formula is the main result of the 
proposed nonextensive generalization of RMT and, to our knowledge, the eigen- 
value distribution functions has never been expressed in this form before. In 
the limit of q = 1, it reduces to eigenvalue distribution functions obtained in 
the RMT and used for further investigation of the Gaussian ensembles. 



3 Nonextensive Wigner surmises 

Analytical expressions for the nearest-neighbor spacing (NNS) distribution can 
be obtained for the 2x2 matrix ensembles. For the Gaussian ensembles, this 
approach leads to the well-known Wigner surmises 

pp(s) = aps 13 exp (-bps 2 ) , (14) 

where (a p ,b fi ) = (f,f), (§, |) and (j|r,f) for GOE, GUE, and GSE, 
respectively. They present accurate approximation to the exact results for the 
case of N — > oo. 

We consider the non-extensive generalization of the Wigner surmises, hoping 
that the results present a reasonable approximation to the physically interesting 
cases of large N, as they do well in the case of Gaussian ensembles. We rewrite 
Eq. (13) for the case of N — 2, and introduce the new variable s = \xi — £2 
and X = (xi + X2)/2 to obtain 

Pp(q,s,X) = C0(q)sP 



l + (q-l)r} q [-s Z + 2X 2 



(15) 



We first consider the case of q > 1 where no limitations are imposed on the 
values of the variables x\ and xi- Thus, we integrate over X from —00 to oo, 
and obtain 

P (q, S) = Ap{q) S P [1 + Bp (q) s 2 ]'^^ , (16) 

where the quantities Ap(q) and Bp(q) can easily be expressed in terms of the 
factors Cp(q) and rj q , and can also be determined from the conditions of nor- 
malization and unit mean spacing as 



Mq) = 7^ / — V l y > ( 17 ) 



r [ -1 - e. + _l 

1,1 2 T q-l 
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and 



Bp{q) 



i 

9-1 



r 



0+1 

2 



1 2 



1 

9-1 



(18) 



For q = 1, Eq. (16) yields the well-known Wigner distributions ^ exp(— as 2 ), 
a = const. As we increase q, the NNS distribution takes the form of a power-law 
type. The quantity Bp(q) diverges when 



= 1 



/3 + 3 



(19) 



This imposes an upper limit on allowed values q. It is equal to 3/2, 7/5 and 9/7 
for the TsOE, TsUE and TsSE, respectively. The variance of this distribution 
is given by 



r (m)r(^)r(-i-| + i i T )r(-2-f + ^ r 



It diverges unless 



r(^)r 



q < <7o 



i 

9-1 



(20) 



1 



/3 + 4 



(21) 



This imposes physical bound on the admissible values of q, because <jp has 
to be finite in order to force condition (5). At values of the entropic index 
exceeding q^, non-extensive statistics does not apply to the random matrix 
model. Figure 1 shows the evolution of the NNS distributions as q increases 
from 1 to goo- The position of the peak of the distribution, which is located at 
ap = y//3/Bp(q) [l/(q — 1) — 1 - /3], moves from Sl = 0.798 to s x = 0.368 for 
TsOE, from s 2 = 0.886 to s 2 = 0.408 for TsUE and from s 4 = 0.940 to s 4 = 
0.671 for TsOE. Neither reaches 0, the peak position of the Poisson distribution 
cxp(— s) of the integrable systems. The figure shows that distributions (16) 
evolve the shape predicted by the corresponding Wigner surmise towards the 
Poisson distribution, but never reach it. The limiting distributions are 



(4+7T 

24s 



? 2 ) 



7 , for TsOE 
for TsUE 



(l+4 s 2)5/2 ! 

1474560 * 4 • for TsSE 



(22) 



(9+64s 2 ) 



TpZ i 



which all tend to s -3 as s — * oo. Recalling the results obtained for the logistic 
map [S] , we may consider these distributions as signaling the edge of chaos where 
the Lyapunov exponents vanish. The Tsallis statistics cannot be extended to 
more regular regimes probably because the degree of mixing then does not justify 
the use of the principle of maximum entropy. 

Now we consider the case of q < 1. For this case, the distribution (9) has to 
be complemented by the auxiliary condition that the quantity inside the square 
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bracket has to be positive. Thus, we integrate over X from —z to z, where 
z - 



1 — i(l — q)rj s 2 . We obtain 



where 



and 



Pp(q, s) = a /3 (q)s f3 [l - bp (q) s 

2[v( 9 )]^ 2 r(2 + i 



a {q) = 



1 

1-9 



r(^)r(f + T ^) 



r (ep) r (2 



1 

1-9 



r(£±i)r(§ + f + T i,) 



(23) 



(24) 



(25) 



Figure 2 demonstrates the behavior of the NNS distributions for the three en- 
sembles as q decreases from 1 to 0. The figure shows that each of the three 
distributions evolves from that of the Wigner surmise towards a shape describ- 
ing a picket-fence spectrum. The limiting distributions are 



Pf)(P,s) 



2 



3^ 



1-^5-*' 



3/2 



for TsOE 



1 



5 2 7 2 tH ' 



3 6 5 2 7 2 



3/2 



3/2 



for TsUE 



for TsSE. 



(26) 



These distributions have their peaks at s = 1 as that of the picket-fence distri- 
bution, but their widths are not small. 



4 Conclusion 

In summary, the NNS distribution obtained from applying Tsallis statistics to 
RMT maps two routes for the transition from chaos to order. One leads towards 
the integrability described by the Poisson statistics by increasing q > 1, but 
ends at the edge of chaos. The other is followed by decreasing q form 1 to 0. It 
leads to a picket-fence type spectrum, such as the one obtained by Berry and 
Tabor TH| for the two-dimensional harmonic oscillator with non-commensurate 
frequencies. 

In order to get a more complete picture of the non-extensive RMT, we need 
to extend our analysis using ensembles of large matrices and including other 
observables as level-number variance X 2 and spectral rigidity A3. Work on 
these lines is in progress. 
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Figure Caption 

Fig.l. The evolution of the NNS distribution for the three symmetry universal- 
ities from the Wigner form towards the Poissonian (both shown by solid 
curves) by varying the entropic index in the range of 1 < q < q x - The 
dashed curves are for distributions having second moments of 2, equal to 
the Poissonian. The dotted curves axe for the values q — q^ explicitly 
given by Eqs. (22). 
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Fig. 2. The evolution of the NNS distribution for the three symmetry universal- 
ities from the Wigner form (solid curves) towards the picket-fence shape 
by varying the entropic index in the range of < q < 1. The dashed 
curves are for distributions having q = 0.5. The dotted curves are for the 
values q = 0, where the distributions are explicitly given by Eqs. (26). 
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